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SUMMARY 


A simplified method for the experimental determination of the 
moments of inertia, product of inertia, and inclination of the princi- 
pal axes; the associated equipment and techniques; and the application 
of this method to a conventional 13 , 000 -pound airplane are described. 
Measurements were made with the landing gear retracted for full and 
empty fuel conditions. The equipment, which consisted primarily of 
knife-edge supports and restraining springs for the pitch and roll axes 
and a single— shaft torsional pendulum for the yaw axis, was designed for 
increased accuracy as well as for simplicity of operation and ease of 
handling as compared with previous methods. At no time was it necessary 
to hoist or jack the airplane in an abnormal fashion. 

Analysis showed the maxi mum possible error of the inertia measure- 
ments to be ±1.7, ±1.2, and ±0.6 percent of the true moments of inertia 
about the X, Y, and Z axes, respectively. For each suspension system, 
measured moments of inertia of known masses of simple form agreed 
within 0.5 percent with the calculated values. 

The results of brief tests have indicated that suitable application 
of the torsional pendulum would 'permit evaluation of the inclination of 
the principal axes to within less than ±0.1°, which corresponds to an 
error of less than ±35 slug- feet squared in the product of inertia of 
the test airplane. 


INTRODUCTION 


The dynamic— stability problems accompanying the unusual config- 
urations and the increases in the relative density of modern aircraft, 
and the application of rational design procedures to servomechanism 
installations necessitate an accurate knowledge of the dynamic response 
characteristics of the airplane. These response characteristics, in 


2 


MCA TN 2201 


turn, are dependent upon the accurate evaluation of product of inertia 
(or inclination of the principal axes) and moments of inertia. The 
effect of product of inertia was usually neglected in earlier dynamic- 
stability work, but in recent years this effect has become more important 
and can no longer be ignored (references 1 and 2). 

The practical problems involved in the experimental determination 
of the moments of inertia and product of inertia have become critical 
with modern aircraft. It has been the practice to suspend and swing the 
airplane as a compound and as a bifilar pendulum, and to correct the 
resulting data for the displacement of the axes of oscillation from the 
body axes through the center of gravity of the airplane (references 3 
to 6). It is difficult to find a structure from which airplanes weighing 
over 10,000 pounds can be suspended for swinging. Usually a building 
with sufficient strength and space to permit swinging has such a high 
overhead structure that the hoisting and handling problems become unrea- 
sonable as accuracy considerations for a compound pendulum require the 
shortest possible pendulum lengths. Even with short pendulum lengths, 
which might be obtained by hoisting and swinging the airplane high above 
the hangar floor, the corrections required for the transfer of axes 
alone (as shown by the data in references 4 to 7) would be as high as 
200 to 700 percent of the final results. Hence, it can be seen that the 
accuracy of the results using such a swinging system would be dependent 
upon small differences in large numbers. 

In view of the structural, hoisting, handling, and accuracy problems 
involved, it appeared impractical to extend these swinging methods to 
the larger and heavier aircraft of the present and future. These prac- 
tical difficulties led the Cornell Aeronautical Laboratory to employ a 
system of pivots and springs to measure the moment of inertia about the 
pitch axis of a B-25J airplane (reference 7). 

When the problem of measuring the moments of inertia and product of 
inertia of a 13 , 000 — pound airplane arose, it was decided to design and 
install equipment that could be used to measure the moments of inertia 
of aircraft weighing up to 20,000 pounds. This equipment was to be so 
designed that the axes of oscillation would be on or as near as possible 
to the body axes of the airplane. The necessity of hoisting and swinging 
the airplane high above the hangar floor was to be eliminated. Handling 
problems were to be reduced to the point where only the handling and 
jacking techniques such as normally used for checking retractable landing 
gear would be employed. The equipment was to be flexible in principle 
to allow its use on any modern aircraft with a minimum amount of special 
fittings . 

A description of this moment— of— inertia gear and its application to 
the measurement of the moments of inertia and product of inertia of a 
13,000-pound airplane are given in this report. During activities not 
associated with this program, the torsional pendulum was damaged prior 
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to the completion of an accurate determination of the airplane product of 
inertia. Rather than delay a flight program scheduled for the test air- 
plane, the inertia measurements were discontinued. However, sufficient 
product— of— inertia data were obtained to warrant discussion at this time. 


SYMBOLS 


Refer to figures 1 and 2 for clarification of the definition of 
certain symbols. The notation of reference 8 was used as a basis for 
the symbols used in this report. 


A aspect ratio of the surface 

C x , Cy static spring constants of the restraining springs 

for the X— and Y— axes oscillations, respectively, 
pounds per foot 



Cz equivalent spring constant of torsional pendulum, support- 

ing roof truss and airplane cradle, foob-pounds per 
radian 

Dp dihedral— angle correction factor 

plan— form tapeiwatio correction factor 

IX, Iy* Iz moments of inertia about the roll, pitch, and yaw axes, 

respectively (the axes are further defined by subscripts 
ref, prin, etc.), slug— feet squared 

IXZ product of inertia, slug— feet squared 

Ix 0 . moment of inertia about an axis in the XZ plane, 

parallel to the axis of oscillation, inclined from the 
X— body reference axis by an angle 9 , and passing 
through the airplane center of gravity, slug— feet 
s quared 


IZ G moment of inertia of the Z— axis torsional— pendulum gear 

(includes the pendulum shaft and the airplane support 
cradle ) 


Lf 

fuselage length, feet 

•* 

Lx 

perpendicular distance 
ing spring to the X 

from center line of the restrain- 
axis of* oscillation, feet 
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perpendicular distance from the center line of the 

restraining spring to the Y axis of oscillation, feet 


period of oscillation, seconds 

area of the surface denoted by the subscript, square feet 
total volume of airplane, cubic feet 
airplane weight, pounds 

span of the surface denoted by the subscript, feet 
local chord, feet 


geometric average depth of the fuselage, feet 

acceleration due to gravity, 32.2 feet per second per 
second 

vertical component of the distance from the X and Y 
axes of oscillation, respectively, to the airplane 
center of gravity, feet 

coefficient of additional mass of an equivalent flat 
rectangular plate 

coefficient of additional moment of inertia of an equiv- 
alent flat rectangular plate 

coefficients of additional mass of an equivalent fuselage 
ellipsoid for motion along the Y and Z axes, respec- 
tively 

coefficients of additional moment of inertia of equivalent 
fuselage ellipsoid about the Y and Z ' axes, respectively 

perpendicular distance in the vertical plane from the 
X axis of oscillation to the centroid of the side area 
of the fuselage, feet 

component of distance in the XY plane of the perpendic- 
ular distance between the Y axis of oscillation and 
the centroid of the top area of the fuselage, feet 


mean aerodynamic chord 



mean chord of the surface 
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If 

r z 

perpendicular distance in the vertical plane from the 
Z axis of oscillation to the centroid of the side area 
of the fuselage, feet 

Hy 

component of distance in the XY plane of the perpendic- 
ular distance between the centroid of the horizontal- 
tail area and the Y axis of oscillation, feet 


perpendicular distance from the Z axis of oscillation 
to the centroid of area of the vertical tail, feet 


perpendicular distance from the X and Y axes of oscil- 
lation, respectively, to the airplane center of gravity, 
feet 

m 

mass * s l' u € s 

w 

geometric average width of the fuselage, feet 

€ 

angle in the XZ plane between the X— body reference axis 
and the X principal axis, positive when the reference 
axis is nose up, degrees 

9 

angle between the X— body reference axis and an inclined 
axis in the XZ plane, positive when the reference axis 
is nose up, degrees 

P 

air density at test altitude, slugs per cubic foot 
Subscripts 

add mass 

additional mass 

fuse 

fuselage 

ht 

horizontal tail 

knife edge 

axis of oscillation 

meas 

as measured (uncorrected for transfer of axes, additional 
mass, etc.) 

prin 

principal axis 

\ 

ref 

body reference axis passing through airplane center of 
gravity 
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vt 


vertical tail 


wing 


wing 


e 


axis in XZ plane inclined from body reference axis by 
an angle 0 


l 


load condition 1 


2 


load condition 2 


APPARATUS, TESTS, AND ANALYSIS METHODS 


Moments of Inertia About Roll and Pitch Axes 


Of the methods of measuring moments of inertia considered, the most 
promising from the practical and the accuracy standpoints appears to be 
a system whereby the airplane is pivoted about an axis of rotation 
located on the airplane structure and restrained from rotating about 
this axis by a spring. The moment of inertia about the axis of rotation 
is then a function of the spring constant, the location of the spring, 
and the period of the resulting oscillation. The apparatus used in these 
tests, the manner in which the tests were carried out, and the method of 
data analysis are described below. 

The position of the airplane center of gravity was determined by 
weighing the airplane in a tail— up and tail— down attitude while holding 
a known reference point on the airplane at a fixed height. By geometry, 
the horizontal and vertical positions of the center of gravity with 
respect to this reference were calculated from the weight and balance 
data. The positions of the airplane center of gravity for the full and 
empty fuel conditions (load conditions 1 and 2, respectively) are given 
in Appendix A. A sketch of the airplane showing the center— of— gravity 
positions and other pertinent dimensions is given in figure 1. 

Roll axis.— The airplane as set up for measuring the moment of 
inertia about the roll axis (X axis) is shown in figure 3. The two knife 
edges fixing the axis of oscillation were located in the plane of 
symmetry below and astride the center of gravity. The restraining springs 
were attached outboard on the front wing spar. As the knife edges were 
below the center of gravity, the springs were preloaded to stabilize the 
airplane in roll. 

The hoisting and jacking of the airplane necessary to position it 
for testing were reduced to a minimum. The airplane was towed into posi- 
tion with the main landing wheels rolling up on low ramps. The restrain- 
ing springs were secured and the tail was raised. Knife edges mounted an 
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hydraulic Jacks were then positioned under the V-block fittings attached 
to the airplane structure (fig. 4). The hydraulic Jacks were used to 
raise the airplane to permit retraction of the landing gear, after which 
the Jacks were lowered, positioning the airplane for test. This pro- 
cedure was reversed to remove the airplane from the test position. 

A standard NACA position recorder coupled with a l/lO— second timer 
was connected to the left wing tip. An oscillation was induced manually 
at the spring and photographic records of the time histories of the 
resulting oscillations were obtained. A double exposure of the oscilla- 
tion in roll is shown in figure 5. 

The moment of inertia about the axis of oscillation is given by 

Ix knife edge = ^ ~ Wh *) (l) 

where P is the period of oscillation, and C x is the sum of the static 
spring constants of the two springs. The moment of inertia at a given 
test attitude about a roll axis through the airplane center of gravity 
and parallel to the axis of oscillation as obtained from the measured 
moment of inertia about the knife-edge axis is given by the equation 

Ix ref = Ix knife edge - Ix add mass - (g + Vp ) ly 

where Ix add la the moment of inertia due to the apparent addi- 

tional mass effect of oscillation in a fluid medium (air) and the term 
[ (w/g)+Vp]l x 2 represents the transfer of axes and the buoyancy and 
entrapped air corrections. Moments of inertia were measured about two 
axes in the plane of symmetry; one parallel to the body reference axis 
(0=0°), and one inclined from the body reference axis (0=7.60°). Use 
of these equations in the evaluation of the product of inertia is dis- 
cussed later. 

Pitch axis.— A knife-edge and restraining— spring method similar in 
principle and handling procedures to that for the roll axis was used to 
measure the moment of inertia about the pitch axis (Y axis). The air- 
plane as set up for test is shown in figure 6. The V blocks were fas- 
tened to the rear wing spar aft of the center-of -gravity position, and 
the restraining spring was secured to the arrester-hook structure at 
the tail. The same knife edges on hydraulic Jacks as used for the roll- 
axis measurements were employed. The knife-edge and V-block assembly is 
shown in detail in figure r J. The same instrumentation that was described 
for the roll-axis measurements was attached to the tail of the airplane 
to obtain photographic time histories of the oscillations. 

The moment of inertia about the pitch axis passing through the 
airplane center of gravity is given by the equation 
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lyraf - ( c y V - “v) (&) ~ ^ mss - (i + To ) V 2 (3) 

Moment of Inertia About Yaw Axis 


A long pendulum length, so detrimental in the case of the compound 
pendulum, has a favorable effect in the case of a bifilar torsional 
pendulum. It can be shown that the accuracy of the bifilar— torsional- 
pendulum method is increased as the ratio of the suspension length to 
the distance between the bifilar supports is increased. It seemed logi- 
cal to extend this principle to the point where a single-shaft torsional 
pendulum would be used for the yaw-axis oscillations; the axis of the 
pendulum shaft then would be the axis of oscillation. The single— shaft 
torsional pendulum has been used extensively in the past to measure the 
moments of inertia of small objects such as projectiles, missiles, and 
dynamic wind-tunnel models. 

The torsional pendulum shown in figures 8 and 9 was maie from a 
4 . 5 — inch— outside-diameter chrome-molybdenum steel tube with solid end 
fittings. The upper fitting was rigidly secured to a suitable roof 
truss. The lower fitting was connected to the airplane support cradle 
by a pin Joint in such a manner that the airplane was free in pitch but 
restrained in roll by the bending of the shaft and in yaw by the twist- 
ing of the shaft. The legs of the cradle were bolted securely to pri- 
mary structure of the airplane. Slots in the cradle beams permitted 
fore-and-aft adjustment of the legs to allow for various center-of- 
gravity positions. The restoring force was provided by the twisting of 
the shaft. The moment of inertia about the axis of the shaft is given 
by the equation 



where C z is the equivalent spring constant of the system. This 
torsional— pendulum spring constant in foot-pounds per radian was eval- 
uated by measuring angular deflections resulting from known applied 
torques . 

It should be noted that equation (4) is rigorous only when the axis 
of oscillation is a principal axis. When the axis of oscillation is not 
a principal axis, there is a coupling between the rolling and yawing 
motions and complex equations relating the two degrees of freedom must 
be considered. Preliminary estimates for the test airplane (verified 
later by the test results) indicated that the principal axis was dis- 
placed less than 4° from the axis of oscillation. Calculations showed 
that the effects of the rolling on the period of the oscillation in yaw 
would be negligible, so that equation (4) was a valid approximation in 
the present tests. 
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The pendulum, was checked with a test frame having a mass and moment 
of inertia about the vertical axis similar to the airplane to be tested. 
Blocks of lead were added to the test frame so as to increase the moment 
of inertia about the vertical axis approximately 30 percent. It was 
found that the moment of inertia of the lead blocks as measured by the 
pendulum agreed within 0.40 percent with the calculated moments of inertia. 
A photograph of the calibration test frame on the pendulum is shown in 
figure 8 . The moment of inertia of the torsional pendulum and cradle 
about the axis of the shaft was determined experimentally. 

The airplane handling procedures were somewhat more complicated than 
those used for the knife-edge measurements. The cradle legs were bolted 
to the airplane and the airplane towed into position under the pendulum. 

The airplane was then lifted in a level attitude to join the legs to the 
cradle beam. The pendulum length was predetermined such that the dis- 
tance the airplane was lifted was just sufficient to permit landing— gear 
retraction. The cradle was then adjusted so that the axis of the shaft 
was coincident with the Z— body reference axis. The airplane as set up 
for oscillating about the yaw axis is shown in figure 9 . 

The photographic recording instruments were attached to the tail of 
the airplane to measure the yawing oscillation. A torque was applied to 
the airplane and held until any undesirable motion had been damped out. 

The torque was then abruptly released and the airplane oscillated about 
the yaw axis. A double exposure showing the motion of the oscillation 
in yaw is shown in figure 10. The moment of inertia about the yaw refer- 
ence axis passing through the airplane center of gravity is given by the 
equation 


I Z 

^ref 


I 7 - I 7 
Tieas ^add mass 



( 5 ) 


Inclination of and Moments of Inertia 
About the Principal Axes 

Product of inertia and inclination of the principal axes .— It is 
assumed that the vertical plane passing through the center line of the 
airplane is a plane of symmetry. Hence the pitch axis is a principal 
axis, since it is perpendicular to the plane of symmetry, and, conse- 
quently, the products of inertia Ixy and- IZY will be zero. In 
figure 2, let X re f and Z re f be the body reference axes, Xg and Z q 
be a set of axes inclined from the body axes by a known angle 9, and 
the axes Xp r ^ n and Zp r p n be the principal axes, inclined at an angle 
e to the body reference axes. Then the moment of inertia about the Xq 
axis is given by the equation 
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IXg = / z* 2 dm= /( z cos 0 — x sin 0) 2 dm 

= /z 2 cos 2 0 dm — 2 /(xz sin 0 cos 0) dm + 
/x 2 sin 2 0 dm 

I Y = Iv cos 2 0 + In sin 2 0 — 

A 0 A ref ^ref 

2 IXZ re f sin 9 cos 0 


so that the product of inertia referred to the body reference axis is 

IX r ef cos2 9 + I Zref sin2 6 ~ J X0 


IXZ r ef 


2 sin 0 cos 0 


( 6 ) 


Since, by definition of principal axes, ^xZ-nrin e ^ ua -*- s zero > from 
figure 2 P 

IXZ prin = /(x"z") dm = /(z cos e - x sin e) (z sin e + 

x cos e) dm = cos e sin e (/z 2 dm — Jx 2 dm) + 

cos 2e /xz dm = 0 


or 


Ixz prln 




sin 2e 


+ I XZ ref (cos 2e) 


= 0 


hence 


or 


tan 2e 


^ ^XZ r ef 

P ^ref Px ref 



^ P XZ r ef 
P ^ref P ^ref 


( 7 ) 


If the moments of inertia Ix re f and. Iz re f adout the bod y reference 
axes and a moment of inertia Iyg about an axis inclined 0° from the 
X body reference axis are measured, then the product of inertia IXZref 
may be determined from equation (6), and the inclination e of the prin- 
cipal axes with respect to the body reference axes can be determined 
from equation (7). 

The noticeable rolling motions which occurred during the torsional 
swingings suggested another method of determining e . This method is 
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based on the fact that application of a pure yawing moment to the system 
produces no rolling when the axis of oscillation corresponds to a pri n- 
cipal axis of the suspended body. The angle between the reference axis 
and the axis of no roll represents the inclination of the principal axis 
of the airplane and gear combination. It can be shown that the correc- 
tion which must be applied to yield e for the airplane alone is closely 
approximated by the expression 

(Ae) G ~ _ 

Zref ^ref 

where IxZq is composed of the product of inertia of the gear about its 
own center of gravity and the terms involved in correcting for the dif- 
ference in centei^of— gravity location of the gear, airplane, and airplane- 
gear combination. 

Brief tests of a preliminary nature were made for load condition 1 
with the airplane suspended with the X reference axis at various angles 
from +3*7° to —2.9° to the horizontal. A position recorder was attached 
to the left wing tip to measure the amplitude and period of the roll and 
the same instrumentation as used for the yaw axis swingings was used to 
measure the corresponding yaw amplitude and period. 

Principal moments of inertia .— The moment of inertia about the 
Y re f axis will be a principal moment of inertia, hence 

lY ref = lY prin < 8) 

From figure 2, 

Iy . = fz" 2 dm = f(z cos e — x sin e) 2 dm 

= f z 2 cos 2 e dm — 2 fx z sin e cos e dm + 

/x 2 sin 2 6 dm 



or 

IXp r in = IX re f cos2 e + I Z re f sin2 e “ 2 IXZ re f sin e cos 6 (9) 


and 


I Zprin = / x " 2 ** =f(z 
= / z 2 sin 2 e dm 


sin e + x cos e ) 2 dm 
+2 f (xz sin e cos e ) dm + 


fx 2 cos 2 e dm 


or 


I Zprin = I X r ef sin2 6 + I Z r ef cos2 e + 2 I XZ re f sin e cos € 


( 10 ) 
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RESULTS AND DISCUSSION 


Moments of Inertia About Body Reference Axes 


Basic data .— The dimensions and physical characteristics of the 
airplane are given in Appendix A. All measurements of length were made 
several times to at least the nearest 0.01 foot. The airplane was 
weighed eight times and the average values were used in the determination 
of the horizontal and vertical positions of the airplane center of 
gravity. It is believed the positions of the center of gravity were 
known within ±0.02 foot. The period data as obtained from the knife- 
edge and pendulum measurements are given in table I. The period values 
for each run are averages of about 30 cycles for the X and Y axes 
and 15 cycles for the Z axis. The timing error was less than 0.01 
second per minute. A mean value of the period of oscillation for each 
set of runs was used for the determination of the moments of inertia. 

The equipment used for the measurements of the moments of inertia was 
tested by oscillating known masses; the calculated and measured values 
agreed within 0.50 percent in all cases. 

Corrections to basic data .— Additional mass and buoyancy effects 
were considered. Additional mass corrections were made according to 
reference 8 and are included in the sample calculations given in 
Appendix B for load condition 1. The resulting true moments of inertia 
about the body reference axes are given in table II. 

Precision .— The effect on the momen1>-of— inertia calculations of 
the possible errors in the various measured and computed quantities is 
summarized in table III, which shows the percentage error in the true 
moments of inertia due to individual errors in each variable taken one 
at a time. The possible errors in the variables were estimated on the 
basis of the present test techniques and the previous experience of 
references 4 , 6, and 8. The total of the individual percentage errors 
is a measure of the over— all precision of the method. The values of 
±1.7, ±1.2, and ±0.60 percent for Ix re f> ^Yref* 811(1 I Z re f^ respec- 
tively, are slightly lower than the values of ±2.5, ±1.3> 811(1 ±0.8 
estimated in reference 6 for the usual swinging methods. Detailed com- 
parison with the data of references 4 and 6 indicates that, in general, 
the errors in measured moments of inertia are slightly greater for the 
new method than for previous methods because of the direct effect of 
errors in the evaluation of the spring constant C. However, this 
disadvantage is more than offset by the reduction, due to the shorter 
suspension lengths, in the magnitude and resultant errors of terms 
involving transfer from the axes of oscillation to axes through the 
center of gravity. This is illustrated by the fact that for the present 
tests the maximum difference between the measured moments of inertia 
and the true moments of inertia about axes through the center of gravity 
is less than 16 percent of the latter, compared with the 200— to 
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700— percent differences inherent with compound— pendulum methods pre- 
viously used (references 4 to j). 

Inclination of the Principal Axes and 
Principal Moments of Inertia 

Inclination of the principal axes tv two— suspension method.— The 
calculations , from equations (6) and (7), of the inclination of the 
principal axis from measurements of the moments of Inertia about the Z 
and X reference axes and an inclined axis in the XZ plane are given 
in Appendix C. With regard to precision, the net effect on Ixz r f °f 
a small error in the directly measured quantity 9 is small, as is the 
effect of an error of +0.6 percent (see table III) in I Z re f However, 
the term Ix re f cos 2 9 — Ix^ in equation (6) represents the small dif- 
ference between large numbers, so that Ixz re f * s very sensitive to 

errors in Ixref and I xg. Since cos 2 9 is nearly equal to 1, the 
possible error in the difference is approximately equal to the error in 
Ix ref — which arises, in turn, from errors in P, L, and l. 

Table III indicates that these items can cause an error of about ±0.55 
percent in each lx value, giving a possible error in Ix re f cos 2 9 — Ix e 
of about 1.10 percent of Ix re f Substitution of this error in Appendix C 
yields a maximum possible error in Ixz re f of ±659 slug-feet squared 

corresponding to about ±1.84° in terms of e. Computations have shown 
that in order to obtain reasonable accuracy in the analysis or prediction 
of the dynamic lateral stability characteristics of high-performance air- 
planes, it is often necessary to know e to less than ±1.0° (reference 2). 
The accuracy of the two— suspension method could be increased somewhat by 
measuring Ix@ at large angles of inclination 9. However, this pro- 
cedure does not appear promising, in view of the handling difficulties 
which might be encountered with airplanes of large size or unusual con- 
figuration. 

Inclination of principal axes by the "Null" method .— In the 
torsional— pendulum swingings with the airplane X reference axis at 
various angles to the horizontal, there was a rolling motion at all test 
attitudes, so that the inclination of the principal axis of the airplane- 
gear combination was not determined directly. However, as shown in fig- 
ure 11, this inclination could be established by interpolation from a 
plot of the value of the dimensionless ratio of maximum rolling— motion 
amplitude to the corresponding yawing-motion amplitude, where these 
amplitudes were measured across the envelope of the oscillations. The 
data indicate an inclination of 2.4°, with a precision of about ±0.1°. 

The torsional pendulum was damaged prior to measurement of the product 
of inertia of the gear itself. However, it was estimated that the 
correction to e due to the gear would be of the order of +O.3 0 , so 
that e for the airplane alone would be about 2.7°. It is believed 
that, with minor modifications to apparatus and technique, e for the 
test airplane could be evaluated to within ±0.1°, which corresponds to 
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an error of ±35 slug— feet squared in Ixz re f* Although the estimated 
value of e of 2.7° is in excellent agreement with the value of 2.77° 
determined by the two— suspension method, even this must be considered 
as fortuitous in view of the possible error of ±1.8k for the latter. 


Principal moments of inertia .- The principal moments of inertia 
were determined from equations (8), (9) > and (10). The sample computa- 
tions (based on the two-suspension Ixz re f <lata) for load condition 1 
are given in Appendix C. The principal moments of inertia, product of 
inertia, and inclination of the principal axes for the two load condi- 
tions are summarized in table IV. Since e is so small, the moments of 
inertia about the principal axes and the resulting possible errors are 
nearly the same as the moments of inertia about the reference axes and 
the corresponding possible errors (tables II and III). 


Comments on Apparatus and Procedures 


Compared with previous methods, the simplicity of the apparatus and 
the handling procedures cannot be stressed too highly. Handling of the 
airplane was reduced to a minimum and at no time was it necessary to 
hoist or jack the airplane in unnatural positions or to any great height. 
In view of the apparent simplicity and accuracy of the Null method for 
determining the inclination of the principal axes, provisions in the 
Z— axis support cradle to facilitate continuous and accurate changes in 
airplane attitude would be desirable. Since the amplitude of small 
rolling motion is of importance in this method, sensitive roll measuring 
instruments based perhaps on strain gages or an optical lever should be 
employed. 

These methods of inertia measurement can be applied, of course, to 
other airplanes, even to very heavy airplanes, if adequate provision is 
made for increasing the weight— carrying capacities of the loaded members. 
There appears to be no great difficulty in the application of the X-axis 
and Y— axis equipment to other airplanes; the detail suspension design 
would be dependent upon the particular airplane configuration and struc- 
ture. It may be necessary, in some cases, to account for the effect on 
the spring constant of the flexibility of the structure between the 
pivots and spring anchors. 

For the Z— axis measurements the application of the overhead tor- 
sional pendulum is limited by the load— carrying capacity of the available 
supporting structure (design load of present equipment was 20,000 pounds). 
This limitation might be overcome by a torsional pendulum which supports 
the airplane from below. Preliminary estimates indicate the practica- 
bility of such a system which would employ a platform flush with the 
ground as the support cradle. The airplane would be supported from this 
platform at the axle axis of the extended landing gear. 
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CONCLUDING REMARKS 


The methods employed in the present investigation for measuring the 
moments of inertia of a 13 , 000 — pound airplane reduced the handl in g 
problems and inherent inaccuracies of previous methods and appear suit- 
able for extension to inertia measurements on very heavy airplanes. 

The test equipment was checked by measuring moments of inertia of 
known masses; the calculated and measured values agreed within 0.50 per— 
cent. Analysis of the precision of the airplane inertia measurements 
showed the maximum possible errors to be ±1.7, ±1.2, and ±0.6 percent of 
the true moments of inertia about the X, Y, and Z axes, respectively. 

At no time was the maximum difference between the measured moments of 
inertia before correcting for additional mass, transfer of axes, etc., 
and the true moments of inertia greater than 16 percent of the true 
moments of inertia, as compared with the 200- to 700 -percent differences 
inherent in the swinging methods previously employed. 

The airplane product of inertia and inclination of the principal 
axes were determined by two methods. The first method was dependent upon 
values of moments of inertia about an inclined axis in the XZ plane and 
about the X and Z reference axes, and was characterized by possible 
errors of ±1.9° in the derived value of e. The other method utilized 
the coupled motion between roll and yaw which occurred when the airplane 
was yawed about an axis other than a principal axis. Brief tests with 
this method indicated that e could be evaluated to within ±0.1°, which 
corresponds to an error of approximately ±35 slug— feet squared in the 
product of inertia of the test airplane. 


Ames Aeronautical Laboratory, 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics, 
April 5, 1950. 
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APPENDIX A.- PHYSICAL CHARACTERISTICS OF AIRPLANE AND INERTIA 
GEAR AS USED FOR MOMENT-OF -INERTIA MEASUREMENTS 


General 

Type: Single-engine, propeller-driven, two-place dive bomber 

Weight and balance 


Load condition 1 
Basic airplane 

Pilot and observer (A00 pounds) 
Research instrumentation 
23 gallons oil 

300 gallons gasoline (fuel tanks full) 


Weight 13,090 lb 

Longitudinal center-of -gravity position 

Gear up . 30.28$ M.A.C. 

Vertical center-of-gravity position from fuselage 

reference (thrust) line — 0.130 ft 

Load condition 2 

(Load condition 1 less fuel) 

Weight 11,525 lh 

Longitudinal center-of-gravity position 

Gear up 27*12$ M.A.C. 


Vertical center-of-gravity position from fuselage 
reference (thrust) line 


Dimensions for inertia measurements 
X axis 

Perpendicular distance from the axis of the 


spring to axis of oscillation, L x 

0 = 7 . 60 ° 


Static spring constant of the restraining 

springs, C x (total) 

Perpendicular distance from the axis of oscilla- 
tion to the airplane center of gravity, 

Load condition 1 

0 = 7 , 60 ° 

0 = 0 6 


Load condition 2 

0 = 7 . 60 ° 

0 = 0 ° 

Vertical component of the distance between the 
X axis of oscillation and the airplane center 


of gravity h x 
Load condition 1 
Load condition 2 


10.30 ft 
10.21 ft 

5832 lb/ft 


1.3^ ft 
1.93 ft 

1.63 ft 
2.19 ft 


1.93 ft 
2.19 ft 
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Y axis 

Perpendicular distance from the axis of spring 

to the axis of oscillation, L 

Perpendicular distance from axis of oscillation 
to airplane center of gravity, 2y 

Load condition 1 

Load condition 2 

Static spring constant, C v 

Z axis 

Equivalent spring constant of the torsional pen- 
dulum and airplane support cradle combina- 
tion, 

Moment of inertia of torsional pendulum and air- 
plane support cradle combination about axis of 
pendulum shaft, 

Vertical component of the distance between the 
Y axis of oscillation and the airplane center 
of gravity, h y 

Load condition 1 

Load condition 2 

Wing 

Area, S 

Span, b 

Aspect ratio, A 

Taper ratio, K 

Mean chord, c 

M.A.C 

Wing volume 

Dihedral angle (top surface front spar) 

Additional momenh-of— inertia coefficient for X 
swinging of a flat rectangular plate (for A=5.87, 

fig. 4, reference 8), k* 

Tapei>-ratio correction factor (fig. 6, reference 8), D^ 
Dihedral correction factor, (fig. 5, reference 8) Dp 
Distance aft from leading edge of wing to leading edge 
M.A.C 

Fuselage 

Fuselage length, Lf 

Geometric average width, w 

Geometric average depth, d 

Fineness ratio of equivalent fuselage ellipsoid .... 

Width-depth ratio, w/d 

Perpendicular distance in the vertical plane from the 
X axis of rotation to the centroid of side area of 
the fuselage. If 

9 = 7.60 .... 

0 = 0 

Component of distance, in the XY principal plane, of 
the perpendicular distance between the Y axis of 
rotation and the centroid of top area of fuselage. If 


16.49 ft 


3.064 ft 
3.411 ft 
5820 lb/ft 


82,000 ft— lb / 

radian 


216 slug— ft 2 


0.751 ft 
1.005 ft 


422 ft 2 
49.72 ft 
5.87 
2.32 
8.48 ft 
109.3 in 
636 ft 3 
6 ° 


0.88 

0.78 

0.80 

0.03 ft 


34 ft 
3-5 ft 
6.11 ft 
7.04 
1.745 


1.59 ft 
2.51 ft 


0.05 ft 
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Perpendicular distance in the vertical plane from 
the Z axis of rotation to the centroid of side 
area of fuselage, Zf z 

load condition 1 


If 

If 


Zl' 

Z 2 ' 


load condition 2 


Coefficient of additional mass of equivalent fuselage 
ellipsoid for motion along the Y and Z axes 
(fig. 7> reference 8), 



Coefficient of additional moment of inertia of equivalent 
fuselage ellipsoid about the Y and Z axes 
(fig. 8, reference 8), 


'f 

*f 


Fuselage volume 


2.44 ft 
2.72 ft 


1.54 

0.57 


0.44 

1.25 

727 ft 3 


Horizontal tail 

Area, S 107.4 ft 2 

Span, b 19.04 ft 

Aspect ratio, A 3*37 

Taper ratio, \ ... . 2.30 

Mean chord, H 5.65 ft 

Volume horizontal tail 40.3 ft 3 

Component of distance in the XY plane of the fuselage 
of the perpendicular distance between the centroid of 
the horizontal— tail area and the Y axis of rota- 
tion, Zty 16.07 ft 

Coefficient of additional mass of an equivalent flat 

rectangular plate of A=3«37 (fig. 3> reference 8) , k O. 87 6 

Vertical tail 

Area, S 45.7 ft 2 

Span, b 7.78 ft 

Aspect ratio, A 1.32 

Taper ratio, X 2.00 

Mean chord, "c 5.88 ft 

Volume vertical tail 17.94 ft 3 

Perpendicular distance from the centroid of area of the 
vertical tail to the Z axis of rotation, 

Z t , load condition 1 19.31 ft 

lt Zz t load condition 2 19.60 ft 

Additional mass coefficient of an equivalent flat rec- 
tangular plate of A=1.32 (fig. 3> reference 8), k . . 


0.65 


X axes 
(a) 


(*>) 


Y axis 


APPENDIX B.- CALCULATIONS OF MOMENTS OF INERTIA ABOUT 
BODY REFERENCE AXES FOR LOAD CONDITION 1 

Additional Mass Corrections 


9=0°, Zf x =2.51 feet 


(ladd mass) 0 =o° = P -j]“g ^k* D^DpS 2 b^ 


r / \ 2 1 

, + ! kf Lf wd i If ] 

wing [ y ^ lx J Jfuse 


= (0.002378) 


TO (0.88)(o.78)(o.8o) (422) 2 (if9.72)+(i.54)(3k)(3.5)(6.n)(2.5i) 2 "! 


= (0.002378) ( 318,244 + 705^.3^) 


(ladd mass)0_QO = 773*56 slug— feet 2 


0=7.60°, If = 1.59 feet 

A 

then 


(ladd mass)0_y S0 o= 763*^7 slug— feet c 


ladd mass = P 


5 k'fy. L f wd 


L f 2 3d 2 

~T" + ~2n 


fuse 


: f z L f wd 


fuse 


r« s 2 

r k — 


It, 


ht 


I — 1 
VO 
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ro 

o 


= 0.002378 


^ (0.44) (34) (3.5) (6.11) 


(34) 


f + 3(6.ll) : 

2it 


[0.57 (34) (3-5) (6.11) (0.05 ) 2 ] + 

5 (0 ' 8T6) %]of (l6 - 07)2 ] } 

0.002378(19631.85 + i.o4 + 107,635.46) 


Z axis 


I add 


mass 


= 302.64 slug— feet 2 


■'■add mass 


.= P 


5 k ’ f z L f wd \~T + W _ 


Lf 2 3v g 


fuse 


.4 t 


= 0.002378 


0O 2 L 


i (i. 25)(34)(3.5)(6 .ii) 
5 


kf y Lf vd ( if, 


(34) 2 3(3.5) : 

4 2n 


[1.54 (34)(3.5)(6.ll)(2.44) 2 ] + 


r (0.65) ^77^ (19.31) 


=0.002378 (53608.45 + 6663.88 + 51099. 48) 


Iadd mass = 264.85 slug-feet 2 


+ 

fuse 
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Momenta of Inertia About Body Reference Axes 
Through the Airplane Center of Gravity 


X axes 

(a) 0=0° 

From equation (2) 


Ex ref Ix meas Iadd mas3 


| + V(p) 

© 


l x z 


-T 

2 V 


- C xV(^) -laid mass (^) 2 -^(P) (^v ) 2 


¥ (V 2 

g x 


^ref 


= 5832 (10. 2i) 2 (o. 03065) - 13090 (i.93)(o. 03065) - 773.56 - 406.52 (1.93) 2 - 
1421 (0.002378)(1.93) 2 

= 18633.72 - 774.33 - 773.56 - 1514.25 - 12.59 = 18633.72 - 3074.73 

= 15559 slug— feet 2 


(b) 0=7-60° 

From equation (2) 


*0 


= - Uh x(^) 2 -I 


add IMSS -5 


= 5832 (10. 30 f (0.02891) - (13090)(1. 93)(0. 02891) - 763.47 - 4o6.52 (1.34) 2 - 


( 1421 ) (0. 002378 )(1. 34 ) fi 


ro 
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= 17887.11 - 730.37 - 783 .^7 - 729-95 - 6.07 = 17887.11 - 2229.86 


ro 

ro 


Y axis 


Z axis 


Ijte = 15657 slug— feet s 


From equation (3) 


^ + V(p) 
2 


<V‘ 


T Y re f “ I Y m eas X add mass 

' °y V (f) -*7 (fij - 3hlM - 5 (l y )2 - T <’> ( V )2 

= 5820 (16.49 ) 2 (0.01906) - (13090)(0. 751X0.01906) - 302.64 - 406.52 (3.064) 2 
1421 (0. 002378X3. 064) 2 

= 30163.88 - 187.37 - 302.64 - 3816.41 - 31.72 = 30163.88 - 4338.14 


l ref 


= 25826 slug-feet" 


From equation (5) 


J ref 


" Iz meas Iadd masB IZr 


= C 


- I. 


= 82000 (0.44502) - 264.85 
I Z ref = 36011 slug— feet 2 


z v 2it 7 ~add mass 
216 = 36492 - 480.85 


- I 


*0 
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APPENDIX C.- CALCULATIONS OF THE INCLINATION OF THE PRINCIPAL AXIS AND 
THE PRINCIPAL MOMENTS OF INERTIA FOR LOAD CONDITION 1 


Product of Inertia 


From equation (6) 


I XZ re f - 



sin 2 0 + l Y cos 2 0 
A ref 

2 sin 0 cos 0 



when 

0=7. 60° 

then 


1X2 = 13.6011) (0.017^9) + (15559) (0.98252) -15657 = 629.8,3 + 15287.03 - 15657 

ref 0.26219 0.26219 


Ixz re f = 991.11 slug-feet 2 


Inclination of the Principal Axis 


From equation (7) 


e 


= tan 1 — ref = i tan 1 

^vef ~ -'-Xref ^ 

= \ (5. 536°) 


2(991. ll ) 
36011 - 15559 


= — tan 1 
2 


1982.22 

20452 



e = 2.768° 


259.86 

0.26219 


O.O9692 


ro 

uo 
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Principal Moments of Inertia 


ro 

•p- 


From equation (8) 

I Yprin = I Y re f = 25826 slug-feet 2 
From equation (9) 

Ix prin = Ix ref COs2 C + Iz ref BlD ^ 6 “ 2 Ixz ref Sln e 008 e = (l 5559 )( 0 . 99846 )‘ f + ( 3601l)(0. 05551 ) 2 - 
2(991. n)(o.o555i)(o. 99846) = 15511 + 111 - 110 

It = 15512 slug— feet 2 

■^-prin 

From equation (10) 

l Zprin = I x ref sin2 e + I Z ref cos2 e + 2 I xz ref 8in 6 cos € = ( 15559 )( 0 . 05551 f + ( 36011 )( 0 . 998 46 ) 2 + 
2(911. n)(o.0555l)(0. 99846) = 48 + 35900 + 110 

I 7 = 36058 slug— feet 2 

■°prin 
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TABLE I.- PERIODS OF OSCILLATION, MOMENT— OF— INERTIA MEASUREMENTS 

[Load condition 1. Fuel tanks full, crew of two (400 pounds), 
airplane weight 13,090 pounds] 



X axis (roll axis) 

Y axis 
(pitch axis) 

Z axis 


0=7.60° 

0 = 0 ° 

(yaw axis) 

Run 

period 
( sec) 

period 

(sec) 

period 

(sec) 

period 

(sec) 

1 

1.0691 

1.1016 

0.8681 

4.1898 

2 

1.0714 

1.1016 

.8676 

4.1928 

3 

1.0690 

1.0994 

.8676 

4.1856 

4 

1.0692 

1.1000 

.8687 

4.1972 

5 

1.0700 

1.1014 

.8686 

4.1848 

6 

1.0683 

1.0999 

.8665 

4.1952 

7 

1.0690 

1.0996 

.8651 

4.1980 

8 

1.0660 

1.0992 

.8661 

4.1925 

9 

1.0682 

1.1009 

.8668 

4.1840 

10 

1.0692 

1.1011 

.8689 

4.1945 

li 

1.0661 

1.1007 

.8675 

— 

12 

1.0690 

1.0973 

.8680 

— 

13 

1.0655 

1.1017 

.8680 

— 

14 

1.0686 

1.1003 

.8677 

— 

15 

1.0682 

1.0976 

.8651 

— 

16 

1.0693 

1.1004 

— 

— 

17 

1.0676 

1.1008 

— 

— 

18 

1.0672 

1.0975 

— 

— 

19 

1.0715 

1.1007 

— 

— 

20 

1.0681 

1.1008 

— 

— 

21 

1.0667 

1.0972 

— 

— 

22 

1.0696 

1.1013 

— 

— 

23 

1.0680 

1.1002 

— 

— 

24 

1.0657 

1.0983 

— 

— 

Mean period, 
second 

1.0684 

1.1000 

.8674 

4.1914 

Maximum varia- 
tion from mean 
period, percent 

.290 

.254 

.2653 

• 1765 

Maximum varia- 
tion from mean 
moment of 
inertia, percent 

.578 

.330 

.613 

• 351 
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TABLE I.- CONCLUDED 

[Load condition 2. Fuel tanks empty, crew of two (400 pounds), 
airplane weight 11,525 pounds ] 



X axis (roll axis) 

Y axis 

Z axis 

Run 

0=7.60° 

0=0° 

(pitch axis) 

(yaw axis) 


period 

(sec) 

period 

(sec) 

period 
( sec) 

period 

(sec) 

1 

1.0399 

1.0585 

0.8669 

4.1064 

2 

1.0392 

1.0569 

.8661 

4.1171 

3 

1.0403 

1.0585 

.8663 

4.1197 

4 

1.0378 

1.0578 

.8663 

4.1140 

5 

i.o4oo 

1.0585 

.8665 

4.1257 

6 

1.0382 

1.0585 

.8651 

4.1199 

7 

i.o4oo 

1.0565 

.8659 

4.1110 

8 

1.0384 

1.0584 

.8654 

4.1145 

9 

i.o4oi 

1.0581 

.8656 

4.1113 

10 

1.0381 

1.0588 

.8650 

4.1210 

ll 

— 

1.0595 

— 

— 

12 

— 

1.0587 

— 

— 

Mean 

period, 

second 

1.0392 

1.0582 

.8659 

4.1161 

Maximum 
variation 
from mean 
period, percent 

.106 

.161 

.116 

.236 

Maximum 
variation from 
mean moment 
of inertia, 
percent 

.185 

• 351 

.267 

.187 
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TABLE II.- MOMENTS OF 
THROUGH AIRPLANE 


INERTIA ABOUT BODY AXES 
CENTER OF GRAVITY 


Item 

Load 

condition 1 

Load 

condition 2 

I Xq> slug-feet 2 
(0=7.60°) 

15,657 

14,687 

Ix ref> slu S- feet2 
' (0=0°) 

15,559 

14,022 

lY ref> slug-feet 2 

25,826 

25,329 

IZref* slug-feet 2 

36,011 

34,710 



/ 
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TABLE III.- RESULTS OF PRECISION ANALYSIS 


Variable 

Possible error, percent of true moment 
of inertia (+ ) 

Symbol 

Estimated error (± ) 

Ix ref 

^ref 

■*-Zref 

C 

0.5 percent 

0.59 

0.58 

0.50 

L 

0.01 foot 

C\J 

• 

.1^ 

— 

P 

0.0005 second 

.11 

.13 

.02 

ladd mass 

10 percent 

• 57 

.13 

.08 

w 

5 pounds 

<.0l 

<.01 

— 

l 

0.02 foot 

.20 

O 

OJ 

• 

— 

vp 

10 percent 

<.01 

.01 

— 

0 

0.01° 

— 

— 

— 

Total 

±1.71 

±1.19 

± .60 
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TABLE IV.- PRINCIPAL MOMENTS OF INERTIA, PRODUCT OF 
INERTIA, AND INCLINATION OF PRINCIPAL AXES 


Item 

Load 

condition 1 

Load 

condition 2 

I Xprin-’ slu S~ feet2 

15,512 

14,215 

I Yprin^ slug-feet 2 

25,826 

25,329 

I Zprin^ slug-feet 2 

36,058 

34,517 

IXZ ref , slug-feet 2 

991 

-1155 

e , angle between 
principal axis and 
X reference axis, 
degrees 

2.77 

-3.19 
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Figure /. - Sketch of test airplane showing pertinent symbols for 
momen t- of- inertia measurements 
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Figure 3 .- Moment of- inertia gear, X axis (roll), 0=0°. 
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(c) Front V— block. 


Figure 4 .— X axis knife-edge fittings. 








Figure 5.- Double exposure showing oscillations in roll. 
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Figure 7«— Knife-edge and V— block assembly. 
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Figure 8.— Torsional- pendulum calibration test frame. 
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Figure 9 .— Moment of inertia gear, Z axis (yaw), 0=0° 
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Figure 10.— Double exposure shoving oscillation in yav. 
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Attitude of airplane X reference axis to 
horizon , deg 


Figure II. - Ratio of maximum roll to yaw as a 
function of airplane attitude as measured 
with the torsional pendulum. Load condition i. 
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